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Abstract 

a : 

Ab initio calculations are used to predict that a superlattice composed of layers of LaTi03 

and LaNi03 alternating along the [001] direction is a S = 1 Mott insulator with large magnetic 
Ch ■ 

moments on the Ni sites, negligible moments on the Ti sites and a charge transfer gap set by 

the energy difference between Ni d and Ti d states, distinct from conventional Mott insulators. 

£> . Correlation effects are enhanced on the Ni sites via filling the oxygen p states and reducing the 

in 

Ni-O-Ni bond angle. Small hole (electron) doping of the superlattice leads to a two-dimensional 

00 " 

single-band situation with holes (electrons) residing on Ni d x i_ y 2 (Ti d xy ) orbital and coupled to 

<^ ■ antiferromagnetically correlated spins in the NIO2 layer. 
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Recent advances in atomic scale oxide synthesis provide new approaches to creating novel 
electronic and magnetic states in artificially designed heterostructures (l|. Pioneering dis- 
coveries include a quasi two-dimensional conducting electron gas at insulating oxide inter- 
46 L thickness dependent or strain-induced metal-insulator transition in ultra thin 




10j and at inserted monolayers 



11 



reconstruction [16|, |17j and pseudogap 
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15|, orbital 



as well as emergent magnetism 
18J. These discoveries mark important milestones 
towards the goal of using the new synthesis capabilities to engineer correlation effects such 
as magnetism and superconductivity. 

In this Letter, we use ab initio electronic structure methods to demonstrate the viability 
of internal charge transfer and heterostructuring as a route to controlling correlation effects 
in complex oxides. The calculations indicate that a superlattice composed of layers of 
LaTiC>3 (a narrow gap 5 = 1/2 antiferro magnet) and LaNiOs (a wide band paramagnetic 
metal) is a 5 = 1 charge transfer Mott insulator, with magnetic moments on the Ni sites, 
negligible magnetism on the Ti sites, and a charge transfer gap controlled by the energy 
difference between the Ni and Ti d levels, which is distinct from conventional 5 = 1 Mott 
insulators (e.g. NiO and L^NiC^). We further show that small hole doping may lead to a 
two-dimensional single-band Fermi surface, with carriers residing on Ni d x 2_ y 2 orbital and 
strongly coupled to correlated spins on Ni d 3z 2_ r 2 orbital, whereas doping with electrons 
also leads to a two-dimensional Fermi surface but with carriers on Ti d xy orbital and weakly 
coupled to antiferromagnetic ordering in the MO2 layer. 

Charge transfer at oxide interfaces has been previously considered, for example, in the 
context of propagation of magnetic order from a manganite to a cuprate material 19|, |20] , 



or titanate 



211 ]. or nickelate 



22 



24l | . Here we show that the charge transfer can be large 



enough to fundamentally change the electronic properties. In the superlattice considered 
here, the lone electron in the conduction bands of TiC>2 layer is completely depopulated and 
the MO2 layer is transformed from a wide band metal to a 5 = 1 Mott insulator. 

Formal valence changes are also known in related materials, for example, A2BB'Oq double 
perovskites, which can be viewed as a [111] direction superlattice 25l|. Experimental indica- 
tions of formal valence changes have been reported for L^MnNiOg 26| and L^TiCoOg 271 ]. 
but these systems seem not to have been the subject of systematic theoretical studies. Here 
we provide a comprehensive theoretical analysis of this phenomena in the context of a dif- 
ferent and more easily synthesized [001] superlattice. An important feature of the [001] 



superlattice not found in double perovskites is an orbital ordering 28|, [29] with Ni d 



r 2_,,2 



at the valence band maximum and Ti d xy at the conduction band minimum. This orbital 
ordering leads to a strongly two-dimensional character for doped holes and electrons which 
are coupled to correlated spins in the MO2 layer. 

To understand the origin of charge transfer and formal valence changes, we sketch the 
band structure of (LaTi03)i/(LaNiOs)i superlattice in Fig.[D\. We begin with bulk LaTiC>3, 
where the nominal electronic configuration involves one electron in Ti d states {d 1 ), which 
lie well above the filled oxygen p states. Bulk LaNiC>3, on the other hand, has a nominal 
electronic configuration of Ni dJ with six electrons filling up the ti g shell and one electron in 
the e g orbitals, which overlap in energy with the oxygen p states. The different electroneg- 
ativities of Ti and Ni suggest that in a superlattice, the Ti d electron may be transferred 
to the NiC>2 plane, leading to a new nominal electronic configuration: Ti d° + Ni d s , with 
strongly enhanced correlations arising from the half filling of the Ni e g shell. 

To quantitatively test the above idea of engineering correlation effects via charge transfer 
and explore the predicted phenomena, we performed density functional plus U calculations 
(DFT+U) including fully structural relaxations (see Supplementary Material for details). 
The unit cell used in the simulations is shown in Fig. [2j We considered ferromagnetic 
(F), stripe antiferromagnetic (S, wavevector (0, |)) and checkerboard antiferromagnetic (G, 
wavevector (|, |)) states and several values of U. The main conclusions are found to be 
independent of the nature of the magnetic order and the value of U (provided that U is 
within a physically reasonable range). The upper left panel of Fig. [TXJ shows the density 
of states (DOS) of the superlattice calculated using the physically reasonable values U-m = 



6 eV and U^\ = 4 eV [30|; the ferromagnetic phase is chosen for clarity of presentation. 
The bands corresponding to the majority spin Ni d states are fully occupied, while those 
corresponding to Ti d states are empty, indicative of a nominal Ti d° + Ni d 8 configuration. 
The superlattice is insulating, with a 0.4 eV energy gap; the highest occupied states are Ni 
d states; the lowest unoccupied states are Ti d states. The energy gap to the Ni unoccupied 
d states is much larger, around 1.5 eV. 

The other panels of Fig. [Tp show the densities of states calculated within the same scheme 
for three reference materials: face-centered NiO (lower left), cubic perovskite LaTiC>3 (upper 
right, note shift in energy axis) and cubic perovskite LaNiOs (lower right). We see that the 
superlattice DOS strongly resembles that of NiO, except with Ti d states added in the middle 
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FIG. 1: A) Schematic band structures of component materials LaTiC-3 and LaNi03. The dashed 
purple lines are the Fermi levels for the two materials. LaTiO"3 shows insulating behavior with a 
small excitation gap set by Ti d-d transitions and a wide energy separation between Ti d states 
and O p states. LaNi03 exhibits metallic behavior with strong mixing between Ni d states and O p 
states. The red arrow highlights the direction of charge transfer in the superlattice. B) Schematic 
band structure of (LaTiO"3)i/(LaNi03)i superlattice. Ti d states are above the Fermi level (dashed 
purple line). Correlation effects split Ni d states into lower and upper Hubbard bands, separated by 
[/ni- C) Densities of states for majority (above axis) and minority (below axis) spins of superlattice 
(upper left) and reference materials NiO (lower left), LaTi03 (upper right; zero of energy is shifted 
so that oxygen bands align with those of LaNiOs) and LaNi03 (lower right). The densities of 

states are obtained using DFT+f/ calculations with £7ni = 6 eV and U^\ = 4 eV. 
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FIG. 2: Side view A) and top view B) of the unit cell of the (LaTi03)i/(LaNi0s)i superlattice 
simulated here. The atom positions are obtained from a representative relaxed structure found in 
our first-principles calculations. La atoms are green, Ti atoms are blue, Ni atoms are black and O 
atoms are red. The oxygen octahedra shaded orange enclose Ni atoms and the oxygen octahedra 
shaded dark blue enclose Ti atoms. The side view A) shows weak tilting of oxygen octahedra. The 
top view B) highlights a strong rotation of oxygen octahedra. 



of the NiO insulating gap. The densities of states of the two constituents LaNiC>3 and LaTiC>3 
(right panels of Fig. HP; note the shift in Fermi levels) are strikingly different. LaNiC>3 is 
metallic, with a large density of non-Ni (in fact oxygen) states in the vicinity of the Fermi 
surface. LaTiC>3 has partially occupied Ti d states. The dramatic difference in densities of 
states between the superlattice (upper left panel of Fig. HP) and the constituent materials 
(right panels of Fig. HP) shows that an electronic reconstruction has taken place (Note 
that in our calculations LaTi03 is metallic because we have assumed cubic symmetry. A 



similar ca 



culation but using the experimental structure for bulk LaTiC>3 produces insulating 



behavior 3l|. This difference is not important here.). With O p states of the two materials 



aligned (see the right panels of Fig. HP), the lowest occupied Ti d states are about 3 eV 
higher than the highest occupied Ni d states, which is the driving force for the electronic 
reconstruction, consistent with the schematics shown in panel A of Fig. HI We mention here 
that the actual charge transfer, which we have computed from the ^-dependence of the full 



DFT charge density (see Supplementary Material), is only ~ 0.2 electron per Ni rather than 
one electron per Ni as inferred from the formal valence changes, because the charge transfer 
in the near Fermi surface states is to a large degree compensated by rehybridization effects 



similar to those discussed in Ref. 321 ] . 



The electronic reconstruction is also revealed by structural distortions. After relaxation 
of the atomic positions, the apical O atom connecting the Ti02 and MO2 layers is found to 
move about 0.05 A towards the Ti atom, so that the apical Ni-0 bond length is 2.01 A and 
the Ti-0 bond is 1.91 A, as expected if charge is transferred from Ti to Ni. The in-plane 
Ti-0 and Ni-0 bond lengths are, on the other hand, much closer to each other and to the 
average of the apical bond lengths: 1.96 A versus 1.94 A. We therefore have an unusual 
situation of a large Jahn- Teller-type distortion about nominally spherical Ni (5=1) and Ti 
(5=0) ions. In addition to the bond lengths, the bond angles are also different from bulk 



values. In the superlattice, the Ni-O-Ni bond angle is 



bund to be 157°, compared to the 



bulk values of 165° (experiment [33|]) and 168° (theory 34j). The decreased Ni-O-Ni bond 
angle reduces the in-plane inter-Ni hopping and helps stabilize the Mott insulating state. 
Both the bond lengths and bond angles have a weak dependence on U; the numbers cited 
here are obtained from calculations at U^i = 6 eV and C/pi = 4 eV, but the 0.1 A difference 
between the apical Ni-0 and Ti-0 bond lengths and a significantly decreased Ni-O-Ni bond 
angle are found at all the U values studied. 

Within the DFT+[7 approximation for a wide range of U^i, we have studied three locally 
stable states: ferromagnetic (F), stripe (5) antiferromagnetic ordering (in-plane wave- vector 
(0, ^)) and checkerboard (G) antiferromagnetic ordering (in-plane wave- vector (|, |)). Panel 
A of Fig. [3] shows the energy differences between these three states. We see that at all values 
of U^i the G-type antiferromagnetism has the lowest energy. As U-m becomes greater than 
1 eV, the energy differences decrease, consistent with the notion that the MO2 planes are 
in a Mott insulating state for which the magnetic energies scale as J ~ t 2 /U. 

Panel B of Fig. [3] shows the fundamental energy gap of the superlattice (which, as seen 
from the upper left panel of Fig. HP, is the energy difference between Ni d states and Ti d 
states). We see that in all cases, the fundamental gap increases linearly with Um and with 
the same slope, indicating that the nature of the magnetic order only affects the onset of 
insulating behavior and not the basic properties of the state, again consistent with the idea 
that a Mott insulator has been created. Further insight into this phenomenon comes from the 
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FIG. 3: Panels A-C: electronic and magnetic properties of (LaTi03)i/(LaNiO"3)i superlattice for 
different magnetic orderings (F denotes ferromagnetic ordering, S denotes 5-type or stripe anti- 
ferromagnetic ordering and G denotes G-type or checkerboard-pattern antiferromagnetic ordering) 
as a function of £%; with U^i = 4 eV. A) Energy difference per Ni between different magnetic 
orderings. B) Minimum excitation gap. C) Absolute value of site projected magnetization of Ni d 
states. D) Metal-insulator boundaries for each magnetic ordering on the {U^i, Uti) phase diagram, 
respectively. F, S and G have the same meaning as in A)-C). 'M' and T denote metallic and 
insulating states. The open symbols denote the upper limits of Um for which metallic phases are 
found, while the solid symbols indicate the lower limits of Um for which insulating phases have 
been found (note that only integer values of U have been studied). 

densities of states shown in Fig. HI We see that changing the nature of the magnetic ordering 
modifies the details of the highest occupied states, but does not affect basic features such as 
the energy splitting between the Ni d main peak (Ni £ 2 g states) and Ti d states (conduction 
band minimum), as highlighted by the maroon arrow in Fig. HJ Panel C of Fig. [3] shows 
the magnetic moment on the Ni d site (computed using the VASP default atomic projector 
for Ni-see Supplementary Material). We see that for U^i > 2 eV, the magnitude of the 
local magnetic moment is essentially independent of the nature of the ordering, is weakly 
[/-dependent and is close to the value expected for the naive atomic 5 = 1 states, once again 
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FIG. 4: Comparison of densities of states for ferromagnetic (F), stripe antiferromagnetic (S) and 
checkerboard (G) magnetic orderings calculated at Um = 6 eV and C/xi = 4 eV. The conduction 
band minimum is aligned (highlighted by the orange line). The energy gap between the conduction 
band minimum and Ni d main peak (denoted by the maroon arrow) is approximately the same for 
different magnetic orderings. 

confirming the Mott nature of the state predicted here. 



To this point, our calculations have used the physically accepted value U^i = 4 eV [30] . We 
have investigated the robustness of our results to the choice of U by sweeping the phase space 
spanned by U-n and U^ for the F, S and G magnetic states. We find the energy sequence 
G < S < F throughout the phase space. Figure [3p shows that although the position of 
the metal-insulator transition boundaries depends on the U values, the essential features are 
[/-independent. The phase boundary for all the magnetic states has a negative slope, which 
is understood as follows: within DFT+U, Uti increases the Ti ti g states because Ti tig shell 
is less than half filled, while U-Ni increases the energy gap between the Ni e g majority and 
minority spins. Therefore, with a given Utu we need a U^ large enough to separate Ni e g 
and Ti ti g states and open a fundamental gap and vice versa. With Ujh 6 eV or larger, Mott 
physics dominates, as the system is rendered insulating for all the magnetic orderings and 
all the values of U-n. 
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FIG. 5: Orbitally resolved Ni e g and Ti ti g densities of states (red: Ni d x 2_ y 2] blue: Ni d 3z 2_ r 2] 
green: Ti d xy ; purple: Ti d xz + d yz ), computed at Urn = 6 eV and U^\ = 4 eV for ferromagnetic 
(F), stripe antiferromagnetic (S) and checkerboard antiferromagnetic (G) orderings. 



The most important experimental test of our calculations is that the superlattice should 
be a quasi two-dimensional 5 — 1 magnetic insulator, with the moments residing on the Ni 
sites and the correlations most likely antiferromagnetic. The position of the apical oxygen 
(nearer to the Ti than the Ni) and the reduced Ni-O-Ni bond angle are other important 
indicators of the predicted electronic reconstruction. Orbital splitting is also a test of our 
calculations. Fig. [5] presents the orbitally resolved densities of states for a representative 
case (Um = 6 eV and U^\ = 4 eV). In the majority spin channel both Ni d x 2_ y 2 and dj >z 2_ T 2 
orbitals are occupied, leading to a spin 5 = 1 configuration, while all Ti ti g orbitals are 
empty, resulting in a spin 5 = configuration. Since both Ni e g orbitals are filled and all 
the three Ti t 2g orbitals are empty, there is no orbital polarization in the sense of Ref . 35|, |36j • 
However, our first-principles calculations show that the valence band maximum is mainly of 
Ni d x 2_ y 2 character and the conduction band minimum is dominantly of Ti d xy character. 
This orbital ordering indicates that small hole doping may lead to an orbitally selective Mott 



metallic state 



37 



with a single orbital (Ni d x 2_ y 2) active at the Fermi surface but coupled 



to antiferromagnetically correlated spins on the Ni d% z 2_ r 2 orbitals. In this case, the carriers' 



motion is strongly affected by magnetism, which is similar to La2_ x Sr x Ni04. However, on 
the other hand, small electron doping may result in a Ti d xy single orbital Fermi surface, 
in which the carriers are only weakly coupled to the antiferromagnetic ordering in the MO2 
layer. This feature is absent in conventional doped Mott insulators. The comparison of 
carriers' strong/weak coupling to correlated spins in one system poses an interesting open 
question and remains to be explored in experiment. Resonant x-ray photoemission may shed 
light on the fundamental differences between these predicted phases. 

In conclusion, we have shown within the DFT+U approximation that in a 
(LaTi03)i/(LaNi03)i superlattice, correlation effects are greatly enhanced on the Ni sites, 
producing a S = 1 Mott insulator with an unusual charge transfer gap set by Ni d and 
Ti d states, distinct from conventional Mott insulators. The superlattice structure leads to 
a single orbital character at band edges and indicates interesting new physics with carrier 
doping, which is absent in double perovskites. We believe that DFT+U approximation pro- 
vides a reliable description of correlation effects for stoichiometric and ordered materials. 
More sophisticated methods (e.g. DFT + dynamical mean field theory) might be needed to 
treat doped systems with mobile carriers. Our findings open a new direction for the control 
of correlation effects via charge transfer and heterostructuring in perovskite oxides. 

We thank Charles Ahn and Sohrab Ismail-Beigi for helpful conversations. This research 
was supported by the Army Research Office under ARO-Ph 56032 and the United States 
Department of Energy under grant DOE-ER-046169. 
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